Quantum state transfer in imperfect artificial spin networks 
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High-fidelity quantum computation and quantum state transfer are possible in short spin chains. 
We exploit a system based on a dispersive qubit-boson interaction to mimic XY coupling. In 
this model, the usually assumed nearest-neighbors coupling is no more valid: all the qubits are 
mutually coupled. We analyze the performances of our model for quantum state transfer showing 
how pre- engineered coupling rates allow for nearly optimal state transfer. We address a setup of 
superconducting qubits coupled to a microstrip cavity in which our analysis may be applied. 

PACS numbers: 03.67.-a, 03.67.Hk, 75.10.Pq, 42.50.-p, 85.25.Dq 



Many protocols for Quantum Information Processing 
(QIP) assume an arbitrary amount of control over the 
evolving system. However, any external influence may 
introduce errors and decoherence. Moreover, having a 
fine single-qubit control over a large register may be 
hard, especially for closely spaced subsystems. This has 
motivated some recent proposals in order to reduce the 
amount of the control over a quantum computer [1| . A 
register can be controlled by global action on clusters of 
qubits or via operations on just an ancillary control qubit 
that is programmable to any desired dynamics 0]. An- 
other possibility is the engineering of an interaction able 
to accomplish a prefixed task. Thus, a process is for- 
mulated according to the program: initialization of the 
register, evolution by the designed interaction and mea- 
surement without further interference on the dynamics. 

Recently, it has been recognized that Heisenberg and 
XY interactions with nearest-neighbor (NN) couplings 
can be used to transfer a quantum state through a net- 
work of qubits virtually without external control 0- 
For XY couplings, linear chains of three NN coupled 
qubits achieve a perfect transfer fidelity and can be 
used as building blocks for a longer communication wire. 
This strategy is useful in those situations where the use 
of a photonic bus (conventionally accepted as a good 
information-carrier) is not easy or convenient. Moreover, 
a three-qubit chain of its XY coupling allows for univer- 
sal quantum computation Very recently, quantum 
cloning in many-qubit networks has been studied [4j . 

The appealing possibilities offered by even simple net- 
works motivate the research for practical systems in 
which XY interactions with possibly tunable couplings 
can be realized. We use the off-resonant interaction of 
a group of qubits with a common bus in order to sim- 
ulate spin networks of tunable XY couplings. In this 
model, the NN restriction is naturally relaxed and we 
study the efficiency of a state transfer process. The in- 
sight we gain through this kind of simulation under con- 



trollable conditions can then be applied to real systems 
of solid-state physics where the amount of control is in 
general smaller. Properly designing the coupling rates, 
we show that the transfer fidelity in chains of different 
lengths can be nearly optimal. Furthermore, in proper 
conditions, once we set the strengths of the inter-qubit 
couplings in a chain of N elements, the quantum state to 
be transferred can be collected at specific times at any of 
the N—k qubits of the chain (k — 0, 1, .., N—2). This may 
represent an advantage with respect to strategies that 
bypass the chain and connects directly the sender qubit 
to the receiver. In the latter, the couplings have to be 
re-designed each time the receiver's position is changed. 

A setup is addressed combining quantum optics 
and superconducting quantum interference devices 
(SQUIDs) [5j that can embody our model. This choice is 
motivated by the advantages of a strong coupling regime, 
the fixed positions of the qubits in the cavity and the 
long life-times achievable in a high-quality factor (high- 
er) cavity. These features have been recently exploited 
for solid state-quantum optics interfaces |(J. 

The system- We consider N > 2 qubits placed inside a 
cavity providing a single boson mode. We assume that, 
via an external potential we can modulate the transi- 
tion energy E q ^ of the qubits. This is possible in many 
QIP devices such as trapped ions, neutral atoms (through 
Stark fields) and SQUID-based systems (via a magnetic 
flux modulating the Josephson energy y|). 

The free Hamiltonian of the qubit system is Y^iLi Hi = 
(1/2) J2f=i E q ,i{&)d-f where of is the z-Pauli operator of 
the i th qubit. Its eigenstates {|±)}i are the basis for the 
i th qubit. The dependence of the qubit energies on <f> 
is explicitly shown. In many cases, the single address- 
ing of the qubits in a register is a difficult task. We 
thus assume that 4> acts collectively on the qubits that 
are closely spaced in the cavity. The cavity field mode of 
its frequency u a , described by the annihilation (creation) 
operator a (cr), is coupled to the qubits. The wavelength 
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A of the field is taken much longer than the dimension 
d of each qubit and their separation so that any depen- 
dence on the position in the cavity is neglected. This 
assumption can be relaxed if necessary. We consider the 
generally valid field-qubit interaction model 

N 

53£ i ,„ = X;n i (at + o)(* i - + *+) (h = i), (i) 
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where we assume the Rabi frequency fii depends on a 
dimensionless parameter rji which can be designed qubit 
by qubit and set once for all. This is in the spirit of pro- 
gramming the system to accomplish a given task without 
interferences to its evolution. We have introduced the op- 
erators af = ((7j r )t = |< Altogether, the dynamics 
of the system is given by H = J2i=i Hi + H a + J2i=i H %a, 
with H a = ui a a'a. Motivated by the recent achievement 
of high-Q cavities (strip-line cavities with Q > 10 4 are 
relevant to this work) and strong-coupling regime in some 
systems that can embody our model |a,bj, we start with- 
out considering dissipation. The main sources of decoher- 
ence in our proposal will be addressed later in this paper. 

We take fij <C u) a , E q ^ under the rotating wave ap- 
proximation. In the interaction picture and with the 
tunable detunings Si = u> a — E qi (§), the Hamiltonian 
is H = J2iLi Qia^tJ~e lSit + h.c. If Si 3> fii, the qubits act 
as a dispersive intra-cavity medium that changes uj a ac- 
cording to u a — ► u a + J2i No real energy-exchange 
is possible, in this case. This permits the elimination of 
the bosonic mode from the dynamics of the qubits. Now 
the evolution is ruled by an effective adiabatic Hamilto- 
nian (where terms oscillating at frequency Si or higher are 
discarded) having the form of a generalized XY model 

- X {SrftJ + a-aj) = £ Xij (of a? + af aj) 

i<j %<j 

(2) 

with Xij = £li£lj/2Sj and the sum runs over the qubits 
having Si = 5j. A state-independent term has been ne- 
glected 8]. Similar adiabatic: interactions were studied 
elsewhere . In particular Biswas and Agarwal cosidered 
a quantum state transfer in a particle-chain. For longer 
chains, however, their proposal turns out to be experi- 
mentally quite demanding [9j. 

Christandl et al. Q considered an XY model with NN 
couplings showing that perfect state transfer is possible 
in chains of 2 and 3 qubits. However, Eq. J2J is about 
long-range, non-NN interactions and is a different model 
whose properties we study, here, with respect to the effi- 
ciency of quantum state transfer processes. 

Engineered state transfer- We denote \i) = 
{®k=£i \~)k) ® f° r a state where only the i th 

qubit is in |+). and eg) denotes tensorial product. 

Because [H e ,J2iLi^i] = 0, H e can be diagonalized in 
each subspace having an assigned number £ of excited 
qubits. In particular, the subspace with £ = (i.e. all 



the qubits in the |— ) state) is spanned by |0). The 
dynamic of an arbitrary state with just one excited qubit 
is instead confined in the subspace with 1 = 1, spanned 
by the orthonormal basis {||}} (i = 1,..,N). With the 
decomposition U £ {t) = e -^,tt = J^" 1 e- ie **|^>(^| 
we calculate the evolved state of the network. Here, 
H e ^ is the restriction of Eq. to the n s — dimensional 
subspace with a given £ and t\ is the eigenvalue 
corresponding to the eigenstate V*)- 

Let the initial state \in) of the network have the qubit 
labelled 1 in a superposition of |— ) 1 and while all 
the other qubits are in |— ), i.e. \in) — (3 |0) + T I D • 
This state evolves to (3 |0) +"fJ2k=o e_l ^'Pfei with 
the projections pki — 1 1) • F° r quantum state 
transfer, we are interested in the transition amplitude 
(N_\ C/i(t) |1) = J2k=o e ~ leit PkiP*kN' whose square mod- 
ulus gives the probability for the process — > |+)jv 
to occur. In general, the N th qubit at a time t is in a 
mixed state p N = Tr 12 .. N -i{U{t) \in)(in\ W{t)). The fi- 
delity of the process is defined as the quantity T(fi, 7) = 
((-I /?* + (+! 1*)pn(/3 |-)+7 |+)). To cancel dependences 
on the initial state, we average T(fi, 7) over the surface 
of the Bloch sphere as T = (1/4tt) / F{fi, j)dT, with dT, 
the surface element. 
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FIG. 1: Coupling configuration for N = 3, 4. If the qubits are 
connected by all equally coupled, the systems are equivalent 
to triangular and squared clusters. If one connection is broken 
by reducing its coupling rate, the system becomes topologi- 
cally equivalent to a two-qubit chain (b) or, generally, to a 
superposition of inequivalent chains (d). 

We start with a simple case that captures the spirit 
of this study. The topology of a network depends cru- 
cially on the configuration of the couplings in Eq. J2J. 
In general, for equal <Vs and Qj's, the coupled qubits 
form a graph of connected vertices. In Fig. ^ we show 
the cases of iV = 3, 4 that result in a triangular and a 
squared cluster respectively (Fig. ^ (a) and (c)). Each 
vertex and solid line represent a qubit and a non-zero 
coupling, respectively. The natural question raised here 
is the performance of the state transfer in this network 
of connections. To analyze it, we explicitly solve the 
problem formulated above, considering just the relevant 
subspace with 1 = 1 but for a generic number of qubits. 
For simplicity, we set Xij = x. We note n s = N and find 
that H e> i admits N — 1 degenerate eigenvalues e\ = — x 
(fc = 0, ..,N - 2) and the eigenvalue ef " Y = (N - l)x 
completes the spectrum. The degenerate eigenspace 
is diagonalizable and an orthogonal basis can be built 
through the Gram-Schmidt algorithm. We get U\{t) — 
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FIG. 2: (a): T vs. the rescaled time r for N = 3. The 
straight line is the bound for classical transfer. We show T for 
the all-equal couplings case (dotted line) and the engineered 
cases with / = 10 (solid line), / = 5 (dashed line), (b): f 
vs. r for the 1 — > 2 transfer process for / = 1.1 (dotted line) 
and f — 5 (solid line). 



e -i(w-i)*t|^v-i\,^,iv-i| with 
The transfer probability 
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K" 1 ) = (i/V^v)E 

is then | (iV|^i(r) |1) | 2 = (2/A 2 ) [1 - cos(AY)], where 
t == is the rescaled interaction time. The resulting av- 
erage fidelity T is plotted against t in Fig.|21(dottcd line). 
The maximum of this function has to be contrasted with 
2/3, the best fidelity achievable for the transfer of a qubit 
state through a classical channel [ljj. It is clear that T 
may be higher than the classical limit. However, trying 
to extend the system, we find that already for N = 4 
T max ~ 2/3 and the quantum channel becomes useless. 

For N = 3, the unwanted coupling X13 does not com- 
promise the state transfer, even if the fidelity is not op- 
timal. However, a strategy to enlarge the range of N 
for which this quantum process is still worthy is desir- 
able. An intuitive approach is to get rid of the redun- 
dant connections. For example, as shown in Figs. ^ (b) 
and (d), by cutting the coupling 1 <-> N, the clusters 
become equivalent to a chain of three qubits (b) and 
to a superposition of elementary chains (three inequiv- 
alent paths connecting 1 to 4 can be found in Fig. 
(d)). The cut may be realized by properly setting ab 
initio the Rabi frequencies Qj's. This is possible by 
an appropriate choice of each r)i defined to be related 
to f2j in Eq. Q). For N — 3, for example, we take 
fli = fi 3 <C to get X12 — x 2 3 X13 . This reduces 
the complexity of the network to a 3-qubit chain. We set 
£12,23 = x and X13 = x/f, with / > 1. The decompo- 
sition of H e> i can be found analytically with the eigen- 
values e° = -x/f, e\ 2 = {x/2f){\ T ^T+W)- They 
correspond to the states \ip\) = (l/>/2)(— |1) + 12} ) an d 
\iPl> 2 ) = M,2(|l) - (e?' 1 /^) |2) + |3)) with the normaliza- 
tions A/i,2 = a;[2x 2 + (e 24 ) 2 ]- 1 / 2 . We find (3| Ui(t) |1) = 
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e i* that reveals the com- 



petition between the different paths the system can follow 
from the sender to the receiver: the path connecting 1 to 
3 via X13 and the one through 1 2) cx A^ 1 |Vi>-A^ _1 \tpl). 
The transfer fidelity .F(0, 1) = |(3|[>i(£)|l)| 2 resulting 



from this interference effect is 



^(0,1) = - + 2A" 2 A" 2 cos(Ae 21 <)+^A; 4 -A; 2 cos(A e l t) 



(3) 

with Ae\ 3 = e\-e[ (i, j = 0, 1, 2). Eq. © reduces to what 
is found by Christandl et al. [2| when / — > 00. For this 
coupling-engineered system, it is ^"(0, 1) ~ 0.973 for / = 
10 and F(0, 1) ~ 0.898 for / = 5 at t ~ 2fn/y/l + 8/ 2 . 
As time goes by, the interferences lead to collapses and 
revivals of the fidelity. The average fidelity T can be 
computed and it is shown in Fig. (a) where it is seen 
that the classical bound can be beaten. 

Increasing the dimension of the network the diago- 
nalization of Eq. J5J) becomes demanding due to the 
dimension of the subspace with £ = 1. However, a 
recurrence law in / and N for the spectrum of H e 
can be found to allow an analytical expression for the 
transfer fidelity regardless of the length of the clus- 
ter. When xin = x/f is the smallest coupling, we 
find that H e has N — 3 degenerate eigenvalues, the 
rest of the spectrum being given by e^~ 3 = —x/f and 

k = 1,2 and sgn(l,2) = ±. The state \N) 
does not belong to the subspace spanned by the 
degenerate eigenstates. The decomposition of U\ 
can thus be effectively performed just considering 
the eigenstates |V>f~ 3 ) = (1/V2) (- |1) + |AQ) and 
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The fidelity keeps the structure of 
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. As / in- 



0, A/jv-2 - * l/v2 and an A-dependent 
value of t when J-(0, 1) ~ 1 can be found. This is shown 
for N = 4, 5, 6 in Fig. where / = 5 has been taken. 
The plot shows the fast oscillations of J-(Q, 1) around an 
average that is sinusoidal with r. Practically, this may 
be a problem as the value of the fidelity depends on the 
ability of stopping the process at precise moments. The 
instantaneous value of T(Q, 1) thus loses significance in 
favour of a mean value F m defined as the semi-distance 
between the maximum and the minimum of the oscilla- 
tions. Because of the decreasing contribution by the beat- 
ing term cx cosfAe^ 2 t], the amplitude A of the fast 
oscillations decreases as N and / grow. This stabilizes 
the fidelity around the slower behavior. A compromise 
between a fast process and the uncertainty in F(0, 1) is 
desirable. From Fig. El (b) (if and ▲) we see that / = 5 
is a suitable choice, allowing a high fidelity (F m > 0.9 for 
N > 4) within r < 60 (up to N — 10) and nearly optimal 
average fidelity T , as shown in Fig. [3] (c). However for 
/ = 1.1 and 4 < N < 10 we find F m > 0.7 (Fig. El (b), 
■) and T > 0.9. It is noticeable that, even for values of / 
that do not optimize the above analysis, we still get very 
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FIG. 3: (a): Fidelity ^(0, 1) vs. r for N = 4, 5,6 and / = 5. 
(b): Amplitude A of the fast oscillations of J-(0, 1) potted 
against N (left vertical axis) for / = 1.1 (♦) and / = 5 {~k)- 
The right vertical axis shows F m vs. N for / = 1.1 (■) and 
/ = 5 (a), (c): Average fidelity T for the values in (a), 
(d): f of the "1 -> 4" (solid line) and "1 -> 3" (dotted line) 
process for N = 4 and / = 1.1. 



good values of T . Furthermore, we stress that, once we 
set the Xij's allowing for the "1 — > N" transfer, the same 
network can be used, in proper conditions, for " 1 — > j" 
processes too (2 < j < N - 1). This is shown in Fig. 
(b) for N = 3 and and in Fig. (d) for N = 4. 

The proposed setup- As a setup combining strong cou- 
pling regime and fixed qubits positions in the cavity, we 
address a cavity-superconducting qubits system || 0- 
An array of mutually coupled Josephson junctions may 
be used as a high-fidelity quantum channel [llj ]. Our 
approach is different as we exploit the advantages of a 
dispersive bus. We take N SQUIDs 5] in a ID super- 
conducting strip-line resonator (Q > 10 4 , uo a — 10 GHz, 
A ~ 1cm d ~ lfim). The whole setup could be 
fabricated via nanolithographic techniques allowing for 
a precise control and calibration of the characteristics of 
the system. The strip-line cavity minimizes the photon- 
losses and protects, to some extent, the qubits from the 
environment. The dephasing due to charge-coupled, low- 
frequency noise is a dangerous source of decoherenc el 13 . 
However, at the degeneracy point of the SQUIDs (Jj. 
each qubit is encoded in the space spanned by the 
equally-charged states \±) i = (V2)(\0) ± |2e» (2e is the 
charge of a Cooper pair) . The environment is not able to 
distinguish between these states which are less sensitive 
to external charge fluctuations 0, 0] . Some estimates 
puts the dephasing rate in the range of 1 usee Q , in 
a resonant single SQUID-cavity systems. On the other 
hand, for ~ 100MHz and S t ~ 1GHz, the Pureed ef- 
fect in this off-resonant setup enhances the characteristic 
life-times up to ~ 50 /isec, suitable for quantum state 



transfer up to N — 10 (/ = 5). We have solved the mas- 
ter equation describing the evolution of networks of up to 
5 qubits, when the decay of cavity field and qubits is in- 
cluded. We found that the resulting dynamics follow the 
trends shown in Fig. [3] (a) without significant deviations. 
In this setup, the qubits energies E q ^ are the Josephson 
coupling energies and the control $ is a proper magnetic 
flux piercing the whole group of SQUIDs. Small relative 
differences between the detunings [upto 1 MHz) do not 
affect the transfer fidelity. In refs. [a,tP]) qubits and cav- 
ity mode are capacitively coupled and the parameters rji 
that can be set properly designing, at the building stage, 
the capacities between the qubits and the the cavity. The 
initialization of the system can be performed if two values 
of $ can be arranged. For $ = $i, we assume that the 
X ij S 0/1" 6 small enough to turn off the mutual couplings. 
At low temperatures (assumed to guarantee the charg- 
ing regime |5j), the off- resonant coupling to the cavity 
makes the probability that a qubit is in |+) . negligible. 
An electrode coupled to qubit 1, then, provides a pulse 
that prepares the desired state. Switching to the proper 
$2 <C $i, suitable for state transfer, the process begins. 
Setting back $ = <I>i, the interaction can be stopped. 

Remarks- A dispersive qubits-boson interaction is able 
to mimic an XY model with tunable couplings. How- 
ever, unwanted non-NN connections appear, complicat- 
ing the dynamics of the network. We have studied quan- 
tum state transfer in this generalized XY model, show- 
ing that the spoiling effect of the redundant connections 
can be bypassed. We have considered a setup combin- 
ing quantum optics and SQUIDs, identifying the param- 
eters discussed in the theoretical model and addressing 
the suitable strategies for quantum state transfer. 

Acknowledgements- We acknowledge support by UK 
EPSRC, Korea Research Foundation (2003-070-C00024) 
and IRCEP. We thank R. Fazio for discussions. 



[1] 



[2] 

[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 

[10] 

[11] 
[12] 
[13] 
[14] 



D. Janzing est al, Phys. Rev. A 67, 04230 (2003); R. 
Raussendorf and H.J. Briegel, Phys. Rev. Lett. 86, 5188 
(2000); S. Benjamin, New J. Phys. 6, 61 (2004). 

S. Bose, Phys. Rev. Lett. 91, 207901 (2003); M. Chris- 
tandl et a/., Phys. Rev. Lett. 92, 187902 (2004). 
M.-H. Yung et al, Quant. Inf. Comp. 4, 174 (2004). 
G. De Chiara et al, quant-ph/0402071 (2004). 
Y. Makhlin et al, Rev. Mod. Phys. 73 357 (2001). 
M. Paternostr o et al, Phys. Rev. B 69 214502 (2004). 
A. Blais et al, cond-mat /0402216 (2004). 
One-qubit terms are usually assumed to be compensated. 
S.-B. Zheng and G.-C. Guo, Phys. Rev. Lett. 85, 2392 
(2000); A. Biswas and G.S. Agarwal, quant-ph/0305101 
M. Horodecki et al, Phys. Rev. A 60, 1888 (1999). 
A. Romito, C. Bruder and R. Fazio, private comm. 

F. Plastina and G. Falci, Phys. Rev. B 67, 224514 (2003). 

E. Paladino et al, Phys. Rev. Lett. 88, 228304 (2002). 

G. Falci et al, Proc. of the Workshop MS+S2004 (2004). 



